Core polarization for the electric quadrupole moment of neutron-rich Aluminum 

isotopes 



OS 

o 
o 

(N 
00 



o 

=5 
S3 



> 

in 
o 

m 

O 

o 



Kenichi Yoshida 

RIKEN Nishina Center for Accelerator-Based Science, Wako, Saitama 351-0198, Japan 

(Dated: February 18, 2009) 

The core polarization effect for the electric quadrupole moment of the neutron-rich 31 Al, 33 Al and 
35 Al isotopes in the vicinity of the island of inversion are investigated by means of the microscopic 
particle-vibration coupling model in which the Skyrme Hartee-Fock-Bogoliubov and quasiparticle- 
random-phase approximation are used to calculate the single-quasiparticle wave functions and the 
excitation modes. It is found that the polarization charge for the proton Ids/ 2 hole state in 33 Al is 
quite sensitive to coupling to the neutrons in the p/-shell associated with the pairing correlations, 
and that the polarization charge in 35 Al becomes larger due to the stronger collectivity of the 
low-lying quadrupole vibrational mode in the neighboring 36 Si nucleus. 

PACS numbers: 21.10.Ky; 21.60.Jz; 27.30.+t 



I. INTRODUCTION 

The nuclear structure far from the /3-stability line has 
been studied very actively with the development of the 
new-generation radioactive-isotope beam techniques to- 
gether with of the microscopic nuclear models applicable 
to drip-line nuclei carried out by the high performance 
computers. 

The ground state properties and the dynamical prop- 
erties such as low-energy excitation modes and giant 
resonances in the medium-mass to the heavier nuclei 
have been successfully described by the self-consistent 
mean-field theory or the nuclear density-functional the- 
ory (DFT) 0. The nuclear DFT has been applied to 
the exotic modes of excitation in unstable nuclei 0, H[ , 
and developed towards description of nuclei in the whole 
chart. Along this line, the self-consistent random-phase 
approximation (RPA) including the pairing correlation 
and the nuclear deformation has been recently developed 
by several groups 0, H, Hi • 

Presently, the small excitation energy of the 2+ state 
and large transition probability B(E2; + — > 2f ) in 32 Mg 
has been discussed in connection with the breaking of 
the spherical magic number N = 20 in neutron-rich 
systems @, B ■ The recent gyromagnetic- factor mea- 
surement of 33 Al at GANIL pj| and the /?-decay study 
of 33 Mg at NSCL [U indicate that 33 Al has a certain 
amount of 2p — 2h intruder configuration. 

The electric quadrupole moment (Q moment) repre- 
senting the deviation from a sphere is directly related 
to the deformation property of the nucleus, and thus its 
investigation for neutron-rich nuclei at around N = 20 
is str ong ly desired both experimentally and theoreti- 
cally [12| , Quite recently, the Q-moment measurement 
of 33 Al has been performed at GANIL [l]|. 

In order to investigate the ground-state Q moments of 
neutron-rich Al isotopes in the vicinity of the "island 
of inversion" [Til ], we carry out the particle- vibration 
coupling (PVC) calculation based on the Skyrme den- 
sity functional, on top of the self-consistent quasiparticle- 
RPA (QRPA). 



The article is organized as follows: In Sec. (TTJ the 
method is explained. In Sec. IIII1 we perform the nu- 
merical calculations and investigate the core polarization 
for the electric quadrupole moments in 31 < 33 < 35 A1. Sec. lIVI 
contains the conclusions. 



II. METHOD 
A. Microscopic particle-vibration coupling model 

The nuclear Hamiltonian of a PVC model [15] on top of 
the Skyrme-Hartree-Fock-Bogoliubov (HFB) and QRPA 
is written as 

H = J2 E M@i + ^B{B X + H couplc . (1) 

i x 

Here Ei is the quasiparticle energy obtained as a self- 
consistent solution of the Skyrme-HFB equation, /3j,j3i 
the quasiparticle creation and annihilation operators. 
The nucleon creation operator ffi(r) is then represented 
using the quasiparticle wave functions as 



(2) 



The phonon energy hw\ is a solution of the QRPA 
equation on top of the Skyrme-HFB, and B\,B\ the 
phonon creation and annihilation operators. We solve the 
Skyrmc-HFB+QRPA equations in the m-scheme. De- 
tails of the calculation scheme is given in Ref . @ . 

Let us now consider the change of the density g(r) 
due to the collective vibrations as g(r) — > g(r) + 6 g(r , t) . 
The nuclear potential C/[g(r)] is accordingly changed as 
U[g(r)] — ► U[g(r) + Sg(r,t)]. To first order in the change 
of the density, the difference of the nuclear potential is 
evaluated to be 



U[g(r) + S Q (r,t)}-U[g(r)] 



dr' S m^Sg(r',t). 



Sg(r') 



(3) 
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Then, the PVC Hamiltonian has a form of 



where 



iJcoupic^ J drdr' 5 -^0y Sg(r', t)ft (r)^(r) . (4) 

We introduce the vacuum defined by the product of 
the HFB vacuum and the QRPA vacuum; 



hi = <0|AQ/3>> 



ft|0} = 0, 



Ba|0) = 0. 



(5) 



The density variation Sg(r, t) can be written in a second 
quantized form using the QRPA modes as 



(6) 



where 6g\(r) is a transition density to the QRPA state 
|A)=fil|0). 

The $1(3] and fa fa parts in ^ {r')^{r) are taken into 
account in the QRPA phonons [ltij . Consequently, the 
PVC Hamiltonian in the leading order reads 

Couple = E/ drdr' ^j^ [8 Qx {r')B{ + ^(t^a] 

The higher order effects can be treated systematically in 
the Nuclear Field Theory 0. 

The coupling interaction in Eq. |(7J) is derived from the 
Skyrme density functional. In the present calculation, 
we approximate the momentum dependent terms in the 
Skyrme interaction by the Landau-Migdal (LM) form. 
This approximation is made only for the construction of 
the PVC Hamiltonian as in Refs. The isoscalar 

(IS) and the isovector (IV) coupling interactions are ex- 
pressed as 



v T =°(r)6gf{r')6(r-r') 
v T=1 {r)5g l x v (r')5{r - t')t z t' z . 



The explicit expressions for v T °(r 
v T=1 {r) = F^/Nq are given in Ref. [2C 



(8) 

F /N and 



B. Description of odd-yl systems 

In order to describe the odd-^4 nuclear systems, we 
diagonalizc the Hamiltonian fl} within the subspace 
0l\O), B\$] |0)}. Then, the resulting state vector is writ- 
ten as 

i^> = E c ^1io)+e^a/3>}- ( 9 ) 

i \j 

The operator for the quadrupole moment can be writ- 
ten as 

Q=(Q) + J2 Qvftfo + E (Q^l + (io) 

ije-nr A 



dr{3z 2 - r 2 )[<pi > i(r)fi,j(' r ) ~ ( P2,i( r ) t P2,j( r )^ 

(11a) 



Qx = (0\[B x ,Q]\0) = j dr{Zz 2 -r 2 )5 Q l{r), (lib) 

and (Q) is the vacuum expectation value. 

The electric Q moment of the eigenstate \cf>) is then 
calculated as 



zQ\4>) =e{(Q)+E 



E C Xj C XkQjk\- 



(12) 



X,jk 



We apply this model to odd-Z nuclei to calculate the 
proton polarization charge of the state \i) for the Q mo- 
ment, which is defined as 



-"pol 



where (i\Q\i) = (Q) 



(WOW 
V (W) 



Parameters 



(13) 



For the mean-field Hamiltonian, we employ the SkM* 
interaction (2lj | in the present numerical applications. We 
use the lattice mesh size Ap = Az = 0.6 fm and a box 
boundary condition at (ft max = 9.9 fm, z max = 9.6 fm). 
The quasiparticle energy cutoff is chosen at -E qPi cut = 60 
MeV and the quasiparticle states up to fl v = 15/2 ± are 
included. The pairing strength parameter is determined 
so as to reproduce the experimental pairing gap for neu- 
trons in 34 Mg (Aqxp.;/ = 1.7 MeV) obtained by the three- 
point formula [H]. The strength t' = -295 MeV fm 3 for 
the mixed-type pairing interaction with the exponent of 
the density-dependence 7=1 leads to the pairing gap 
(A„) = 1.71 MeV in 34 Mg [|. On top of the Skyrme- 
HFB, we solve the QRPA equation within the space of 
the two-quasiparticle excitation of E a + Ep < 60 MeV. 
The momentum-dependent terms in the residual interac- 
tion are exactly treated. 

The density-dependent Landau parameters Fq and Fq 
for the PVC interaction are determined by the parame- 
ters of the SkM* interaction. As shown in Refs. @, HH , 
the attraction of the LM interaction is stronger than that 
of the self-consistent interaction in which the momentum- 
dependent terms are treated exactly. Therefore, we mul- 
tiply the overall factor /lm for the PVC interaction. We 
use /lm = 0.8 and 0.9 for comparison. Accordingly, the 
difference of the results can be considered as a theoretical 
uncertainty 



III. RESULTS AND DISCUSSION 

A. Properties of 32 ' 34 ' 36 Si 

We describe the odd-Z neutron-rich Al isotopes as a 
proton single-hole state coupled to the neighboring Si 
isotopes because the pairing gaps of protons in Si isotopes 
are zero. We summarize here the ground state properties 
and the structure of quadrupole excitations in 32 > 34 > 36 Si. 

In Table HI the ground state properties are summa- 
rized. The neutron-rich Si isotopes under investigation 
are spherical although the calculated deformation param- 
eters are not exactly zero ((3 2 = 0.02 in 34 Si). This is due 
to the artificial breaking of the spherical symmetry as- 
sociated with the finite mesh size and the rectangular 
box, and thus it is considered as a numerical error. The 
average pairing gaps of neutrons are finite, while those 
of protons are zero. This indicates that the 34 Si has a 
neutron 2p — 2h configuration in its ground state. The 
neutron occupation number of the I/7/2 orbital is 0.31, 



0.78 and 2.21 in 32 ' 34 . 36 Si and that of the 2p 3/2 orbital is 
0.10 in 36 Si. 

Figure [T] shows the response functions for the isoscalar 
(IS), isovector (IV) and proton quadrupole excitations. 
We can see a prominent peak at around 3 MeV in all 
of the isotopes under investigation. The isoscalar transi- 
tion strengths are J B(IS2;0+ -> 2^) = 626,637 and 1137 
fm 4 in 32 > 34 < 36 Si, corresponding to 104, 98 and 162 in 
Weisskopf unit. Besides the low- lying collective 2 + state, 
we can see the giant quadrupole resonances (GQR) at 
around 20 McV and 30 MeV for the IS and IV excita- 
tions. 

The 2f state in 32 Si is mainly generated by the sd-shcll 
configurations of neutron and proton. The microscopic 
structure of the A 71 " = + component is given by the 
two-quasiparticle excitations of (j/2s 1 / 2 ® 1^3/2) with a 
weight of 0.04, {vld 3 / 2 ) 2 with 0.61, and (ir2si/ 2 8> 1^5/2) 
with 0.30. The microscopic structure of the K v = 1 + 
and 2 + components is the same within the numerical 
accuracy as of the A' 71 " = + component because of the 



TABLE I: Ground state properties of 32 > 34 . 36 Si obtained by 
the deformed HFB calculation with the SkM* interaction and 
the mixed-type pairing interaction. Chemical potentials, av- 
erage pairing gaps, root-mean-square radii for neutrons and 
protons are listed. The average pairing gaps of protons are 
zero in these isotopes. The average pairing gap is defined 
(A), = - J drhg/ J dr§. 
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FIG. 1: (Color online) Response functions for the isoscalar 
(IS), isovector (IV) and proton quadrupole excitations in 
32,34,36 The transition strengths are smeared by a 
Lorentzian function with a width of V = 1 MeV. 



spherical symmetry. The strength in the energy region 
15 < hw < 25 MeV exhausts 79.4% of the IS energy- 
weighted sum rule (EWSR) value. The IS strength in 
the low energy region up to 10 MeV exhausts 1.2% of 
EWSR. The IV strength is distributed in a wider energy 
range 20 < Tilo < 40 MeV. The summed strength in this 
energy region exhausts 75.2% of the IV-EWSR value. 

In 34 Si the neutron excitations into the pj-shell become 
appreciable for the 2f state. The microscopic structure 
of the A x = + component is given by the (^ld 3 / 2 ) 2 
excitation with a weight of 0.21, (fl/7/2) 2 with 0.10, and 
(tt2si/2 55 1^5/2) with 0.63. The strength in the energy 
region 15 < hu; < 25 MeV exhausts 79.8% of the IS- 
EWSR value, and the strength in the energy region 20 < 
hw < 40 MeV exhausts 74.3% of the IV-EWSR value. 

The = 0+ component of the 2^ state in 36 Si 
is mainly generated by the (fl/7/2) excitation with a 
weight of 0.39, {v2p 3/2 ) 2 with 0.08, and (tt2s 1/2 ® ld b/2 ) 
with 0.47. The strength in the energy region 15 < hoj < 
25 McV exhausts 77.6% of the IS-EWSR value, and the 
strength in the energy region 20 < fhu < 40 MeV ex- 



hausts 72.9% of the IV-EWSR value. Due to the mixing 
of the IS and IV modes in neutron-rich nuclei, we can 
see an appreciable IV strength in the lower energy region 
10 < hoj < 20 MeV. The summed strength in this energy 
region exhausts 10.4% of the IV-EWSR value. 



B. Polarization charges in 



3 A1 



For describing the r = 5/2+ state of 31 > 33 > 35 Al, we 
diagonalize the Hamiltonian fTJ) in the model space of 
the proton single- hole state of the orbital \{W = 

-5/2+)- 1 ) and the coupled states of \{W = -5/2+)- 1 ® 
wjr=o), K-3/2+)- 1 8>wjr=i) and \{-l/2+)- 1 ® lo k=2 ). 
We take the QRPA states \u>\) whose IS or IV quadrupole 
transition strengths possessing greater than in 1 W.u. 

The dimension of the Hamiltonian ((TJ is 225 for 31 Al. 
The -P = 5/2 + state of | 31 Al) is constructed mainly by 
the one hole state and the hole coupled to the 2+ state 
as 

| 31 Al; ^ = 5/2+ M = 5/2) = 0.93K-5/2+)" 1 ) 

+ 0.21|(-5/2 + ) -1 ® 2+(K = 0)) 



0.23|(— 3/2+) -1 
0.16K-1/2+)- 1 



2+(K = l)) 
• 2+(K = 2)). 



(14) 



The amplitude associated with the other components 
are smaller than 0.1. The coupled states in Eq. (fT4"|) 
correspond to the K^l^^) -1 <8> 2 + ) state in the j- 
schemc representation. The ratios of the amplitudes 

are identical to those of the Clcbsch-Gordan coefficients 
,55 ,55, ,53 ,55, ,51 ,55, 

( 22 2 ° l 22 ) < ( 22 21| 22 ) ^ ( 22 22| 22 ) ' The ^ ad ™P ole 
moment of 31 Al is then calculated using Eq. (p~2|) as 
13.6(13.9)efm 2 , where we use f LM = 0.8 (0.9). From this 
value, the polarization charge of Eq. (fT3|) is calculated as 
ep ol = 1.03e(1.08e). The result of the calculation overes- 
timates slightly the experimental value of the Q moment 
11.2±3.2efm 2 [T^. 

For 33 Al, the dimension of the Hamiltonian ((T|) is 250. 
As in 31 Al, the wave function of 33 Al is written by mainly 
the one hole state with a weight (the squared ampli- 
tude) of 0.90 and the hole coupled to the 2 + state with 
a weight of 0.08. The calculated values of the Q mo- 
ment and the polarization charge are 13.0(13.4)efm 2 and 
e poi = 0.89e(0.96e), respectively. The quadrupole mo- 
ments and polarization charges of 31 A1 and 33 Al are not 
very different. 

The dimension of the Hamiltonian (fT]) is 285 for 35 Al. 
As in the case of 31 Al and 33 Al, the wave function of 
35 Al is written by the one hole state with a weight of 
0.85 and the hole coupled to the 2 + state with 0.14. The 
contribution of the coupling to the 2 + state is larger than 
in 31,33 Al. The calculated values of the Q moment and 
the polarization charge are 14.7(15.1)efm 2 and e^ ol = 
1.12e(1.18e), respectively. 
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FIG. 2: (Color online) Same as Fig. [Q but for the RPA 
strengths in 34 Si. 



In order to see separately the effects of coupling to 
the low-lying modes and to the giant resonances on the 
polarization charge, we diagonalize the Hamiltonian ((T|) 
in the model space containing only the RPA modes with 
their energies larger than 10 MeV. The obtained polariza- 
tion charges for the lrf 5 / 2 orbital are e^ ol = 0.20e(0.22e), 
0.21e(0.24e) and 0.18e(0.20e) in 31 A1, 33 Al and 35 Al, re- 
spectively. These values are close to the systematic value 
obtained in Ref. [H (0.21e, 0.18e and 0.15e in 3 1> 33 > 3 5A1), 
where the microscopic PVC calculations were performed 
including only the giant resonances on top of the self- 
consistent HF+RPA in light neutron-rich nuclei. 

The enhancement of the polarization charge in 35 Al 
is thus due to the strong collectivity of the low-lying 
quadrupole vibrational mode in the core nucleus 36 Si be- 
cause the effects of coupling to the giant resonances are 
not sensitive to the neutron number in Al isotopes under 
investigation. 



C. Effects of the pairing correlations in 33 Al 

We investigate the effects of coupling to the p/-shcll 
in 33 Al. The dominant correlation in the present case is 
the pairing correlation because the core nucleus 34 Si is 
calculated to be spherical at the HFB level. 

Figure [5] shows the IS, IV and proton quadrupole tran- 
sition strengths in 34 Si obtained by solving the Skyrmc 
HF+RPA equations without pairing correlations. The 
collectivity of low-lying states are tremendously weak- 
ened, while the structure of GQR is not very different to 
that obtained by solving the Skyrme HFB+QRPA equa- 
tions shown in Fig. [T] 

The 2+ state is constructed dominantly by the 
(ir2si/2 ® 1^5/2) excitation with a weight of 0.98. The 
strength B(IS2;0+ -> 2+) has only 237 fm 4 . Using the 
RPA transition densities, we diagonalize the Hamiltonian 
((T|) with a dimension of 241. The calculated wave func- 
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FIG. 3: (Color online) Same as Fig. [T] but in 32 Mg. 



tion is mainly generated by the one hole state with a 
weight of 0.97 and the hole coupled to the state with 
a weight of 0.02. This state is thus dominantly described 
by the proton s<i-shell configurations because the 2f state 
is generated by the proton excitation to the 2si/ 2 orbital. 

The resulting Q moment is 11.0(11.4)efm 2 , and the 
polarization charge is e* ol = 0.59e(0.65e). 

The Q moment and the polarization charge in 33 Al are 
quite sensitive to the neutron pairing correlation at TV = 
20. This is similar to the enhancement mechanism of the 
B(E2;0 + ->■ 2f) in 32 Mg because the neutron pairing 
correlation is indispensable for the strong collectivity of 
the 2f state in 32 Mg [25|. 



D. 



! Mg as a core 



the proton quasiparticle of the ld$/2 orbital, and the cou- 
pled states of a sd-shell quasi-proton to the quadrupole 
modes in 32 Mg. In the present calculation with a dimen- 
sion of 276, | 33 A1) is constructed by the quasi-proton of 
the ld 5 /2 level with a weight of 0.85, |7rld 5 / 2 ® %t) with 
0.10, |7r2s 1/2 ®2j f ) with 0.04, and |7rld 3/2 (g)2f ) with 0.01. 
The electric Q moment of 33 Al is then calculated as 12.4 
(12.6) efm 2 . This is consistent with the calculation in 
Sec. ni.B. 

IV. CONCLUSION 

The polarization charges for the electric quadrupole 
moment of the neutron-rich Al isotopes at around N = 
20 have been investigated by carrying out the micro- 
scopic particle-vibration coupling calculation in which 
the coordinate-space Skyrme-Hartree-Fock-Bogoliubov 
and quasiparticle-random-phasc approximation are em- 
ployed to calculate the singlc-quasiparticlc wave func- 
tions and the transition densities. 

It has been found that the neutron pairing correlations 
are crucial to generate collectivity of the 2+ state in 34 Si, 
and that the polarization charge of the proton hole state 
of the ld$/2 orbital becomes small in the absence of the 
pairing correlation in 33 Al. The effect of the neutron 
pairing correlation at TV = 20 on the enhancement of the 
polarization charge in 33 Al is very similar to the enhance- 
ment mechanism of the B(E2) in 32 Mg [2o| . 

Effects of coupling to the giant resonances on the po- 
larization charge are not very different in a small region 
of isotopes, and the low-lying collective modes have much 
effect on the polarization charge. Therefore, the polar- 
ization charge in 35 Al is larger than in 31 > 33 A1 as a con- 
sequence of the stronger collectivity of the 2+ state in 
36 Si. 



The wave function | 33 A1) can also be constructed by 
the quasi-proton coupled to 32 Mg. The pairing gaps of 
neutrons and protons are 1.80 and 1.42 McV. Figure [3] 
shows the response functions for the IS, IV and proton 
quadrupole excitations. At 0.9 MeV, we can see a promi- 
nent peak possessing 2414 fm 4 and 430 e 2 fm 4 for the 
IS quadrupole strength and B(E2), respectively. This 
result agrees well the experimental value @. The 2f 
state is constructed by the two-quasiparticle excitations 
of (vld 3/2 ) 2 with a weight of 0.10, (^l/ 7 / 2 ) 2 with 0.09, 
(7rld 5 / 2 <8>2si/ 2 ) with 0.11, and (7rld 5/2 ) 2 with 0.61. 

We diagonalize the Hamiltonian (fT]) within the space of 
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